J/^jj + Xcj Production at the B Factories under the Principle of Maximum Conformality 
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Under the conventional scale setting, the renormalization scale uncertainty usually constitutes 
a systematic error for a fixed-order perturbative QCD estimation. While, the recently suggested 
principle of maximum conformality (PMC) provides a principle to eliminate such scale ambigu- 
ity in a step-by-step way. Using the PMC, all non-conformal terms in the perturbative expan- 
sion series are summed into the running coupling, and one obtains a unique, scale-fixed, scheme- 
independent prediction at any finite order. As an application, we make a detailed PMC analysis 
for both the polarized and the unpolarized cross sections for the charmonium production processes, 
e + + e~ — > J/ip(ifj') + Xcj with (J = 0, 1, 2), up to next-to-leading order (NLO). Then, we compare 
our predictions with the measurements at the B factories, BABAR and Belle, and the theoretical 
estimations obtained in the literature. Because the non-conformal terms are different for various 
polarized and unpolarized cross sections, the PMC scales of these cross sections are different in prin- 
ciple. It is found that all the PMC scales are almost independent of the initial choice of renormal- 
ization scale. Then, the scale uncertainty for both the polarized and the unpolarized cross sections 
are greatly suppressed even at the present NLO level. It is found that the charmonium production 
is dominated by the J = channel. After PMC scale setting, we obtain a(J/tp + Xco) = I2.I4I2 43 
fb and a(jp' + Xco) = 5.141*'q2 fb, where the squared average errors are caused by the bound state 
parameters as m c , | i?.//^ (0) | and |it^, eJ (0)|. 
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PACS numbers: 13.66.Bc, 14.40.Pq, 12.38.Bx, 13.; 



The double charmonium production processes, e + + 
e~ — > J/tp(ip') + Xcj with (J = 0, 1, 2), play an important 
role for understanding the production and hadronization 
properties of the heavy quarkonium. In addition to the 
non-relativistic QCD (NRQCD) factorization pj, such 
kind of double charmonium production also provides a 

I 



good opportunity to learn the perturbative QCD (pQCD) 
color-singlet mechanisms. 

The processes e + + e~ — > J/ip(tp') + Xcj have been 
studied within the pQCD factorization up to next-to- 
leading order (NLO), cf. Refs.@-0|. Its polarized cross 
section up to NLO can be written as 
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where Ai and A2 stand for the helicities of the outgoing J/ip and Xcj, and the tree-level coefficients 
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Here (Ir stands for the renormalization scale, the charm 
quark electric charge e c = 2/3 and the color factor CV= 
4/3. The dimensionless parameter r — Am^/s, and |P| = 



XHs,Mj /i> ,Ml J )/(2^, 



where 



X(x, y, z) — x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. 

The parameters |i£j/^(0)| and |i?^ cJ (0)| are the radial 
wavefunction at the origin and the first derivative of the 
radial wavefunction at the origin of the bound states J/ ip 
and Xcj, respectively. The expressions for the LO coeffi- 

(/Xfl) can be 
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and the NLO coefficients Bi ^ 
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found in Ref.[6|], for convenience, we put them with some 
reformations in the appendix. 

For calculating all the polarized and the unpolarized 
cross sections, one needs to introduce the renormaliza- 
tion scale fiR. As a conventional choice adopted in the 
literature, the value of [1r is fixed to be an initial value 
Mii'*' which is usually taken to be the typical momen- 
tum flow of the process, e.g. 2m c or \fs/2. However, it 
has been observed that under such conventional scale set- 
ting, one usually obtains explicit renormalization scheme- 
and scale- dependent estimation 2 6] . And the scale un- 
certainty usually provides a large systematic error under 
the conventional scale setting, and even a non-convergent 
pQCD series for some high energy processes. 

Recently, it has been pointed out that the principle 
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of maximum conformality (PMC) provides a possible so- 
lution for eliminating such renormalization scale uncer- 
tainty The main idea of PMC lies in that the 
effective or the optimal renormalization scale is gener- 
ally different from its initial value, which is formed by 
absorbing all the non-conformal terms into the coupling 
constant. One can choose any arbitrary value to be /x™*, 
but the optimal PMC scale and the resulting finite-order 
PMC prediction are both to high accuracy independent 
of such arbitrariness, consistent with the renormalization 
group invariance. The PMC satisfies all self-consistency 
conditions for setting the renormalization scale . Af- 
ter proper procedures, all non-conformal {/3i}-terms in 
the perturbative expansion are summed into the running 
coupling so that the remaining terms in the perturba- 
tive series are identical to that of a conformal theory. 
The QCD predictions from PMC are then independent of 
renormalization scheme, inversely, such scheme indepen- 
dence can be adopted to derive commensurate scale rela- 
tions among different observables After PMC scale 
setting, the divergent "renormalon" series (n!/3"a™) does 
not appear in the conformal series and the convergence of 



the pQCD series can be greatly improved. For example, 
such merits have been shown by a next-to-next-to-leading 
order PMC analysis of the top-quark pair production at 
the hadronic colliders HI G3- As the main purpose of 
the present paper, we shall show that even at the NLO 
level, the QCD prediction after PMC scale setting can 
also have such good features, thus our understanding of 
e + + e~ — > J/ip + Xcj can be greatly improved. 

To set the PMC scales for the double charmonium pro- 
duction process e + + e~ — > J/ip + Xcj, we need to decom- 
pose the NLO coefficient ^(^r^) into Ar dependent 
part (non-conformal part) and /^-independent part (con- 
formal part), i.e. 

= bS^ja. + ft'w'), (3) 

where the non-conformal B^^ 2 (//g lt ) and the confor- 
mal B^ c ^\[i'x t ) can be found in the Appendix. Then, 
following the standard procedure of PMC, the polarized 
cross section (fTJ) can be written as, 



a 1 - I ' n-(n rs "' \ J i ■ {^l^V , I 



with the PMC scale 

^M) = /#* exp (Mk nit )) • (5) 

Because the non-conformal terms -Bj[j^ 2 (/•*#'*) are differ- 
ent for J — 0, f , 2, the PMC scales of these cross sections 
are also different. 



Similarly, the NLO coefficients i?/^™*) of the unpo- 
larized/total cross sections can also be decomposed into 
the non-conformal part B^\fJL R at ) and the conformal 
part Bt icon) (fj} R lit ). After PMC scale setting, the total 
cross section 



al = Ilo {l + (i^Vif)A, + ^""W*)) } 



changes to 

J _ J, PMC\ I J -I , as ^Rft,J ) ^J(con) / init\ 

t - a t \^R-,t,J I Ilo \ H B t [fi R ) 



with the PMC scale 



(G) 



^y=^ nit exp(-B t J(/?) (^ nit ))- (7) 



It is noted that the initial scale dependent logarithmic 



terms in B J £™> (^) and B J t (con} (^ it ) should be ab- 
sorbed into the coupling constant simultaneously with 
the non-conformal terms, so after applying PMC scale 
setting, SjJ™ a 



?J(coTl) and B J t {con) are still at the scale ffi. 



As for the functions for the total cross sections ([6]) , we 
have 
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t 2 Ilo = 0q,o Ilo + 2cx 2 i0 |lo + 2ctq,i |i 

/ R 2(con) 2 I , 0R 2(con) 2 i , 0R 2(ccm) 2 i , 0R 2(con) 2 
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for J = 2. Here, for convenience, we have omitted the initial scale (a*™*) dependence to the coefficients B'^ co "' 
the PMC scale (jUr^j) dependence in the LO ' ; 



(8) 
(9) 

(10) 

(11) 
(12) 

(13) 

(14) 
(A5) 

(16) 
and 



cross sections. 



The summed up non-conformal functions B^ (m™*) 
are usually different from their polarized terms 



(Vr*), then the PMC scales for the total cross sec- 



tions are generally different from the PMC scales for the 
polarized cross sections. Furthermore, because the non- 
conformal terms B'^\^ lt ) are different for J = 0, 1, 2, 
the PMC scales of these total cross sections are also dif- 
ferent. Thus for the considered processes, we have in- 
troduced several PMC scales for the polarized and un- 
polarized cross-sections, which is reasonable and is due 
to their different non-conformal terms. It shall be found 
later that all these PMC scales are almost independent 
of the initial choice of renormalization scale, and then 
the renormalization scale can be greatly suppressed. As 
a consistent cross-check of PMC, it is helpful to show 
whether after PMC scale setting, by summing up all the 
polarized cross section, one can still get the same result 
as that of the directly calculated unpolarized cross sec- 
tion. This, inversely, also can be used as a criteria of 
whether the complex higher-order calculations are right 
or not. 

To do numerical calculation, we set the collision energy 
y/s = 10.6 GeV, the charm quark mass m c = 1.5 GeV 
and ct(\fs) = 1/130.9 (ltij . Using the two-loop a s running 



0.1184 hj, we obtain A 



(71/ =3) 



0.386 



with a(Mz) — u.noi i±±ij, vvc uuiaw ii qcd 

GeV, A^= 4) = 0.332 GeV, and A^= 5) = 0.231 GeV. 
Based on the experimental values for the leptonic width 
of J/ip and ip{2S) and the width of Xc2 to two photons, we 
can inversely determine the radial wavefunction |-R rls (0)| 
at the origin and the first derivative of the radial wave- 





Conventional scale setting 


PMC scale setting 




2 m c 


V^/2 




2 m c 


v^/2 




0? (fb) 


9.23 


6.81 


5.21 


12.14 


12.14 


12.14 


o\ (fb) 


1.01 


0.85 


0.70 


0.99 


0.99 


0.99 


0? (fb) 


1.53 


1.26 


1.03 


1.56 


1.56 


1.56 



TABLE I. Initial scale uncertainty for the total cross sections 
of e + + — > J ji\> + XcJ, where three typical initial renormal- 
ization scales are adopted. It shows that the cross sections 
after PMC are almost independent of ^t^ 1 ', and the residual 
scale dependence due to unknown higher-order {/3i}-terms are 
highly suppressed. 



function at the origin \R n p (0)| through the formulas: 



and 

lx(n?>)->77 " 



4a 2 
9m? 



64a 2 
45m;! 



16a s (2m c ) 
3^ 



16a s (2TO e ) 
3tt 



\Rns(0)\ 2 (17) 



\R np (0)\ 2 - (18) 



More explicitly, u sing the central experimental values for 
the decay widths [17| , we obtain: 

\Rj/4>{Q)\ 2 — 0-855 GeV 3 , (19) 
l-^V (0)| 2 — 0.362 GeV 3 , (20) 



|JJ Xcj (0)| 2 = 0.055 GeV 5 



(21) 



We present the total cross sections before and after 
PMC scale setting in Table |TJ Under the conventional 
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scale setting, = the total cross sections are quite 
sensitive to the choice of renormalization scale. For ex- 
ample, by varying /ig xt from 2m c to yfs/2-, total cross 
sections for the cases of Xco, Xci and Xc2 production are 
changed by about 26%, 16% and 8% respectively. Then, 
the scale uncertainty always constitutes an important er- 
ror source for the conventional estimation. On the other 
hand, after PMC scale setting, we observe that the total 
cross sections remain almost unchanged by varying the 
initial renormalization scale /i™' = 2m c , y/s/2, y^s. 
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FIG. 1. Total cross section versus the initial renormalization 
scale /jg* for e + + e~ -> J/ip + Xcj with (J = 0, 1,2), be- 
fore and after PMC scale setting. Here the short notation 
Conv. stands for the cross section derived from the conven- 
tional scale setting. 



More explicitly, we show the total cross section versus 
the initial renormalization scale /x™* in Fig. (JT|) . It shows 
that, after PMC scale setting, the total cross sections 
are almost flat versus the initial renormalization scale, 
then the renormalization scale ambiguity is eliminated 
even at the NLO level. After PMC scale setting, we 
have resummed the non-conformal {/?o}-terms known at 
the NLO level into the coupling constant, while the LO 
cross sections will be significantly increased, the NLO 
corrections are suppressed to a certain degree. Then, 
the perturbative convergence can be improved, in prin- 
ciple. A comparison of our estimations with the exper- 
imental results are shown in Table |TT] and IIII1 where as 
a comparison, we also put the results of Refs.[H, [g] . Our 
estimation for the total cross section of J/ip + \ c o pro- 
duction is consistent with the experimental result, but 
the i/j' + Xco production cross section is still smaller than 
the central value of the data. Since the data is still with 
large error, a future more accurate measurement shall be 
helpful to clarify this puzzle. As a final remark, there is 
residual scale dependence due to unknown higher order 
{/?i}-terms, which however is highly suppressed 

It is also interesting to show how PMC scale setting 
affects the polarized cross sections. We take e + + — >• 



J/ip + Xcj as an explicit example. For each polarized 
processes, the coefficients for the conformal or the non- 
conformal terms are different due to the enhancement 
or cancellation among different polarizations, then it is 
reasonable that the effective momentum flows of them are 
different. So, in principle, we should introduce different 
PMC scales for each polarized process. 

More explicitly, we put the PMC scales for the polar- 
ized cross section versus the initial renormalization scale 
M ™t i n Fig.©, where the PMC scales for J = 0,1,2 
are presented. All PMC scales for the polarized cross 
sections are almost independent of the initial renormal- 
ization scale, which are the same as the case of the total 
cross sections. Then, the scale uncertainty are also elim- 
inated for the polarized cross sections at the NLO level. 
This also shows that our present calculation for both the 
polarized and unpolarized cross sections are consistent 
with each other. Fig.© shows 

• For e + + e~ — > J/ip + XcO process, we need to intro- 
duce two different PMC scales for (Ai,A 2 ) = (0,0) 
and (1,0) respectively. The contribution from the 
helicity channel (1,0) is about 2 times larger than 
that of the (0, 0) channel. 

• For e + + e~ — > J/ip + Xci process, we need to in- 
troduce three PMC scales. The PMC scales for 
(Ai, A 2 ) = (0, 1) and (1, 1) are equal to each other 

due to the same non-conformal terms -Bqi = 
B\^\ so we only have two different PMC scales. 
The helicity channel (0,1) provides the dominant 
contribution for this process. 

• For e + + er — > J/ip + Xc2 process, we need to intro- 
duce five PMC scales corresponding to five helicity 
channels. The PMC scales for (Ai, A2) = (0, 1) and 
(1,1) are the same due to the same non-conformal 



terms _B, 



203) 
0,1 



B\^\ while other PMC scales are 
different, so we have four different PMC scales. The 
contribution from the helicity channels (0, 0) and 
(1,0) are dominant over others. 

Because of the exponential suppression to the PMC 
scale as shown by Eq.©, the LO cross section for each 
polarized cross section shall be increased to a certain de- 
gree, which are shown clearly in Tables IIV1 Ivl and I VII 
Defining a parameter C, being the ratio of the cross sec- 
tion at NLO level to that at LO level, we find that the pa- 
rameter C are 0.73 and 0.82 for (Ai,A 2 ) = (0,0) and (1,0) 
channels of e + + e~ — » J/ip + Xco under the conventional 
scale setting; which however changes down to 0.08 and 
0.49 after PMC scale setting. For e + + e~ J/ip+Xd,2, 
the conditions are similar. This shows that the pQCD se- 
ries become more convergent after PMC scale setting. 

As a summary, we have applied the PMC to study the 
polarized and the unpolarized cross sections of e + + e~ — > 
J/tp + Xc.j{J = 0,1,2) up to NLO level. After PMC 
scale setting, we find that the final results are scale- and 
scheme- independent and the pQCD series becomes more 
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Belle 
a x 6 >2 (o) [18] 


BaBar 
a x B >2 [19] 


our prediction 


a(J/ip + X co) 
a(J/tp + Xd) 
cr{J/tj) + Xci) 

<j{J/^ + Xcl) + <T{J/i> + Xc2) 


6.4 ± 1.7 ± 1.0 

< 5.3 at 90% CL 


10.3 ±2.51^ 


19 ■)> +2.45+1-05 
iz ' i *-2.19-1.05 
n QQ+0-07+0.09 
u - 33 -0.09-0.08 
1 t- R +0.40+0.09 
1 - JO -0.33-0.09 
9 rc+0.47+0.18 
^• JJ -Q.42-0.17 


cr(ip' + Xco) 
<r(i>' + Xci) 

<T(4>' +Xc2) 
<K^' + Xcl)+^' + Xc2) 


12.5 ±3.8 ±3.1 
< 8.6 at 90% CL 




r 1^+1.04+0.44 
J - ± -0.93-0.43 
n ,,9+0.03+0.04 
u -^-0.04-0.04 
fl RR+0.34+0.04 
U.UU_ Q .22-0.04 

1 r >S+°- 37 +0- 08 
l.UO_ 26-0.08 



TABLE II. Comparison our theoretical predictions for J/ip(tj)') + Xcj {J = 0,1,2) with the experimental results at the B 
factories (in unit: fb). The first error of our prediction is for m c £ [1.40 GeV, 1.60 GeV] and the second error is from the 
uncertainty of experimental results for Tj,r n3 \_+ e + e - and F x ( np )_ >77 . 



FIG. 2. PMC scales versus the initial renormalization scale /Jj? 1 * for the polarized cross sections of the e + + e — > J/if> ± Xci 
processes with J = 0, 1, 2. 





J/^±Xc0 


i> ±Xco 


Belle a x B x "° [> 2] [20] 


16 ±5 ±4 


17±8±7 


Belle a x B x -"[> 2(0)] [18] 


6.4 ± 1.7 ± 1.0 


12.5 ±3.8 ±3.1 


BaBar a x B x <=° [> 2] [191 


10.3 ± 2.5 




Wang, Ma and Chao [5] 


9.5 


4.1 


Dong, Feng and Jia [6J 


8.62 


4.98 


Our result 


12 14+ 267 

1^-11-2.43 


r 14+1.13 



TABLE III. Total cross sections (in unit: fb) for J/tp(ip )+XcO 
production at the B factories. Our estimation together with 
those of Refs. 0, [f| are presented as a comparison. 





conventional scale setting 


PMC scale setting 




LO 


NLO 


sum 


LO 


NLO 


sum 


<o (fb) 


0.869 


0.637 


1.506 


2.891 


0.240 


3.131 


<?% (fb) 


1.454 


1.197 


2.651 


3.025 


1.471 


4.496 



TABLE IV. Polarized cross sections for e + + e — > J/tp + Xco 
before and after PMC scale setting, /i 1 ^ 1 * = y/s/2. 



convergent for both the polarized and the unpolarized 
cross sections. So we gain a more accurate pQCD esti- 
mation than those in the literature. More explicitly, 





conventional scale setting 


PMC scale setting 




LO 


NLO 


sum 


LO 


NLO 


sum 


*i.o (fb) 


0.001 


-0.007 


-0.006 


0.002 


-0.023 


-0.021 


4i (fb) 


0.291 


0.113 


0.404 


0.702 


-0.185 


0.517 


(&) 


0.026 


-0.001 


0.025 


0.063 


-0.058 


0.005 



TABLE V. Polarized cross sections for e + + e — > J/ip + Xci 
before and after PMC scale setting, /i 1 ^ 1 * = y/s/2. 





conventional scale setting 


PMC scale setting 




LO 


NLO 


sum 


LO 


NLO 


sum 


4o (fb) 


0.340 


-0.050 


0.290 


1.132 


-1.723 


-0.591 


<?lo (fb) 


0.210 


0.189 


0.399 


0.438 


0.262 


0.700 


4i (fb) 


0.051 


0.003 


0.054 


0.122 


-0.093 


0.029 


(fb) 


0.026 


0.004 


0.030 


0.063 


-0.041 


0.022 


cr?, a (fb) 


0.002 


0.0003 


0.002 


0.003 


-0.001 


0.002 



TABLE VI. Polarized cross sections for e + ±e —¥ J/tl> ±Xc2 
before and after PMC scale setting, /i 1 ^ 1 * = \fsj2. 



• After PMC scale setting, we obtain total cross sec- 
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tioins 

a(J/i> + X co) = 12.14+2;g fb, 
cr(J/^ + Xci) = 0.99t°;^ fb, 
cr(J/^ + Xc2) = 1.56t°;Sfb 

and 

t(V' / +Xco) = 5.14±i:S| fb, 
<7(V' + Xci) = 0.42lH|fb, 
<7(V' + Xc2) = 0.66l°;gfb, 

where the errors are caused by the input bound 
state parameters as the c quark mass, the wave 
function at the origin or the derivative of the wave- 
function at the zero, and etc. Our estimation for 
the total cross section of J/ip + Xco production is 
consistent with the experimental result, but the 
ip' + Xco production cross section is still far less 
than the central value of the data. Since the data 
is still with large error, a future more accurate mea- 
surement shall be helpful to clarify this puzzle. 

• A comparison of Tables HH HVl IVl and [VI] shows 
that e + + e~ — > J/ip(ip') + Xcj processes are dom- 
inated by J = channel, which are measured at 
the B factories. In contrast, the contribution of 
the e + + e~ — > J/i/}(ip') + Xd,ci channels for this 
process seem to be rather modest even after the 
PMC scale setting. The future super B factory, 
with much higher luminosity, may eventually ob- 
serve these polarization patterns. 

• Since at the present NLO level, not all the pQCD 
convergence of the polarized channels are improved 
after PMC scale setting. Maybe, we really need a 
next-to-next-to-leading order (NNLO) calculation 
to improve the pQCD convergence and make the 
pQCD estimation more reliable. This is similar 
to the case of top-pair production at the hadronic 
colliders, where we have found that after PMC 
scale setting, not only the scale dependence can be 
eliminated but also the pQCD convergence can be 
greatly improved by a NNLO calculation 0, S LT3 ■ 

• As a cross check, by taking all the same input pa- 
rameters, we exactly obtain the results listed in 
Ref.@. It is found that after PMC scale setting, 
even though we have introduced different PMC 
scales for different type of cross sections caused by 
the different non-conformal terms, by summing up 
all the polarized cross sections, we can get the same 
result as that of the directly calculated unpolarized 
cross section for the case of J = and J = 1. This 
shows that our PMC procedure is self-consistent. 
However, there is slight difference for the case of 
J = 2. This shows that the expression for the coef- 
ficient i?Q^Q°"^ listed in Ref.Q may have some slight 
typos. 



• It is found that after PMC scale setting, for the 
polarized cross sections a\ (negative even before 
PMC scale setting) and Cq , we obtain the "incor- 
rect" negative cross section. Similar to the QCD 
pomeron case, this could be the problem of the 
unphysical MS-scheme itself [2l|. It is found that 
under the MS-scheme there is also the same "in- 
correct" sign for the NLO hi ghe st eigenvalue of the 
BFKL equation [22j ■ In Ref. [2l|j , it has been shows 
that the reliability of QCD predictions for the in- 
tercept of the BFKL pomeron at NLO when eval- 
uated using the physical schemes, such as the mo- 
mentum space subtraction (MOM) scheme [23[ or 
the T-scheme based on T — > ggg decay, can be sig- 
nificantly improved 0- 

The "incorrect" sign for a\ and Uq under MS 
scheme may show that those unknown/missed non- 
conformal terms from the NNLO or higher orders 
may have sizable contributions to the LO PMC 
scale and can not be safely neglected. In fact, the 
large K factor for these two polarized cross-sections 
shows their pQCD convergence are terrible. How- 
ever, if transforming to the physical schemes, such 
as the MOM scheme or the T-scheme, our initial 
numerical results show that such "incorrect" sign 
disappears as the same as the pomeron case. This 
indicates an interesting feature of PMC scale set- 
ting, i.e. there is is not only the optimal scale but 
also the optimal scheme. A detailed discussion on 
this point shall be presented elsewhere (23 |. 

• As a final remark: e + + e~ — > J/^(tp') + XcJ pro- 
cesses provide another good example for eliminat- 
ing the renormalization scale ambiguity by using 
the PMC. Thus, after applying the PMC scale set- 
ting, it will not only increase the precision of QCD 
tests, but it will also increase the sensitivity of 
the collider experiments to new physics beyond the 
standard model. 
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After PMC scale-setting, the final results shall be independent 
the choice of the renormalization scale due to the commensurate 
scale relation, however, the conformal terms at each perturbative 
order could be different in different schemes, and then pQCD 
convergence could be different. 



7 



APPENDIX: COEFFICIENTS FOR THE LO AND 
NLO TERMS 



The LO coefficients c J x% \ 2 (r) are 



c i,o( r 

c io( r 
c o,i( r 
c i,o( r 



-1 - fOr + 12r 2 , 
-9 + 14r, 
V6(2-7r), 
— Vor, 

2a/6(1 - 3r), 
V2(-l + 2r + i2r 2 ), 
\/6(-l + 5r), 
V2(ffr-3), 
2%/6(f - 3r), 
-2\/3. 



the relations among the same polarization combination 
(•^1,^2) with different J, the non-conformal NLO coeffi- 
cients Si i and the conformal NLO coefficients 5^ c °") 

Ai,A2 pil^ A 2 

at the renormalization scale hr can be written as q 



^0,0 


- ^0,0 


R 0(/3) 
-°1,0 


R l(/3) 


R l(/3) 
-"0,1 


R 2(/3) 


R l(/3) 
-#1,1 


R 2(/9) 
- -#1,1 


-#1,2 





111 



/4 



2 In 2 I , 



(22) 



2(/3) 
1,0 



^( ln^ + ^-+21n2k23) 



17 



i-fein^ 

12 V s 

T2 1 61n — 



9 

121n2 + 13 
121n2 + 13 



- + 21n2 



(24) 
(25) 
(26) 



Using the asymptotic expressions derived in Ref . [6] and and 



0(con 
0,0 

Q(co7i 
1,0 

l(con 
1,0 

l(con 
-*0,1 

■j^l(con 

B, 

B, 



14 

2 (con 
0,0 

2(con 



J 0,l 



rJ 2(con 
-#1,0 



^1,1 



D 2(con 

-#1,2 



2 1 

•-(4-ln2)lnr- -(46 + tt 2 - 40 In 2 + 33 In 2 2), 



^ In 2 r - ^(139- 104 In 2) In r- ^ fl61 + ^- - ^ln2 + 100 In 2 2 



27 



5 In 2 r + (7 - 2 In 2) In r - 19 + 2tt 2 + 75 In 2 - 21 In 2 2 , 



In 2 r - (46 - 99 In 2) In r - - (616 + 74tt 2 - 1696 In 2 + 303 In 2 2) 



6r 
1 

12 
1 

12 



2 1 

--(4-ln2)lnr- -(64 



101n 2 r + 2(l + 171n2)lnr- -(266 + 7tt 2 - 128 In 2 + 147 In 2 2) 



1041n2 + 33 In 2), 



' / 1 - In 2 r - (22 - 43 In 2) lnr - - (284 + 30tt 2 - 380 In 2 + 159 In 2 2) 



= -- 21n 2 r + -(5 + 8 In 2) lnr - — (291- 8tt 2 + 171 In 2 + 312 In 2 2) 
3 V 6 18 

= — [41n 2 r- (46- 621n2)lnr- ^(274 + 27tt 2 - 3161n2 + 91n 2 2) ) , 



12 



--I 2 In r 



-(1 + 13 In 2) lnr + -(-7ir 2 + 140 - 104 In 2 + 237 In 2 2) 
3 9 



(27) 
(28) 

(29) 
(30) 
(31) 
(32) 
(33) 
(34) 
(35) 
(36) 



2 Recently, the authors have presented their erratum for Ref. 0, 



but we think it still has some typos. We have corrected some 
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typos in the non-conformal {/3o}-terms that are listed in Ref.@| 



by using the explicit relations among the same polarization com- 
bination (Ai, A2) with different J. 



